SCATTERING THEORY FOR JACOBI OPERATORS WITH 
QUASI-PERIODIC BACKGROUND 
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Abstract. We develop direct and inverse scattering theory for Jacobi opera- 
tors which are short range perturbations of quasi-periodic finite-gap operators. 
We show existence of transformation operators, investigate their properties, 
derive the corresponding Gel'fand-Levitan-Marchenko equation, and find min- 
imal scattering data which determine the perturbed operator uniquely. 



1. Introduction 

Classical scattering theory deals with the reconstruction of a given Jacobi oper- 
ator 

(1.1) Hu(n) = a(n)u(n + 1) + a(n — l)u(n — 1) + b(n)u(n), 

which is a short range perturbation of the free one Hq associated with the coefficients 
a(n) = |, b{n) = 0. This case has been first developed on an informal level by 
Case in a series of papers jS] - QU]. The first rigorous results were established by 
Guseinov ^Hjj who gave necessary and sufficient conditions for the scattering data 
to determine H uniquely under the assumption 

(1.2) £|n|(|a(n)-i| + |6(n)|) <oo. 

Further extensions were made by Guseinov [T^l, [2D|, and Teschl |27|. Additional 
details and further references can be found, e.g., in |28|. 

In addition to being of interest on its own, scattering theory can also be used 
to solve the initial value problem for the Toda equation via the inverse scattering 
transform. This has been formally developed by Flaschka Jl] (see also [25] and 
for the case of rapidly decaying sequences) who also worked out the inverse 
procedure in the reflection-less case. Further results and an extension of the method 
to the entire Toda hierarchy were given by Teschl in jJB] and |27j . 

The next interesting problem is to replace the free Hamiltonian H by one with a 
periodic potential. First results in the case of Sturm-Liouville operators have been 
obtained by Firsova in a series of papers (see J3j). For further results, including 
potentials with different spatial asymptotics, and additional references see Gesztesy 
et al. 16 . In the discrete case, the investigation has only recently been started by 
Boutet de Monvel and Egorova [2] and by Volberg and Yuditskii who treat 
the case where H has a homogeneous spectrum and is of Szegd class exhaustively 
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from an operator point of view. Applications to the Toda lattice can be found in 
Bazargan and Egorova and Boutet de Monvel and Egorova 0- 
Finally, let us give a brief overview of the paper: 

Section |3 collects some well-known facts from Riemann surfaces and introduces 
the necessary notation. Section [3] introduces the Baker- Akhiezer function and in- 
vestigates the quasi- momentum map. In the periodic case, where the integrals can 
be explicitly computed, this was first done in |24| . In addition, we characterize 
the second solution at the band edges. In Section 0] we prove existence of Jost 
solutions and use them to characterize the spectrum of the perturbed operator. In 
the periodic case, existence of Jost solutions was first shown by Geronimo and Van 
Assche ^3 and the fact that there are only finitely many eigenvalues in each gap 
was first proven in Cojuhari and later rediscovered in Teschl [251 - Sectional 
introduces the transformation operator and proves the crucial decay estimate on 
its coefficients. This was first done by Boutet de Monvel and Egorova [2] in the 
periodic case under the additional assumption that all spectral gaps are open. We 
fix a problem in the original proof and at the same time simplify and stream- 
line the argument. Section El investigates the scattering matrix. Our main result 
here is the reconstruction of the transmission coefficient from the reflection coeffi- 
cient, which was not known previously, even in the periodic case. Sectional derives 
the Gel'fand-Levitan-Marchenko equation and proves positivity of the Gel'fand- 
Levitan-Marchenko operator. In addition, we formulate necessary conditions for 
the scattering data to uniquely determine our Jacobi operator. Our final Section|Hl 
shows that our necessary conditions for the scattering data are also sufficient. It 
should be mentioned that, due to the lack of continuity with respect to the spacial 
variable n, a significant change in the strategy of the original proof in the continuous 
case from [32] is needed. 

Our approach uses heavily the fact that the Baker- Akhiezer function is a mero- 
morphic function on the Riemann surface associated with the problem. This strat- 
egy gives a more streamlined treatment and more elegant proofs even in the special 
cases which were previously known. In this respect it is important to emphasize 
that, in contradistinction to the constant background case, the upper sheet of our 
Riemann surface is not simply connected and in particular not isomorphic to the 
unit disc. 

2. Quasi-periodic finite-gap operators and Riemann surfaces 

To set the stage let M be the Riemann surface associated with the following 
function 

(2.1) RlH^z), R 2g+2 (z) = [] ( z - E i)> E o < Ei < ■ ■ ■ < E 2g+1 , 

3=0 

g e N. M is a compact, hyperelliptic Riemann surface of genus g. We will choose 
i?2g+ 2 ( 2; ) as the fixed branch 

23+1 

(2.2) <+ 2 W = - II V^E~ 
where is the standard root with branch cut along (— oo,0). 
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A point on M is denoted by p — (z, ±R 2g+2 (zj) = (z,±),z€C,orp = oo±, and 
the projection onto CU {00} by n(p) = z. The points {{Ej,0), < j < 2g + 1} C M 
are called branch points and the sets 

(2.3) n± = {{z,±R\' g %{z)) | z e C\ (j[E 2j ,E 2j+1 }} c M 

3=0 

are called upper, lower sheet, respectively. 

Let {a,j, bj}j =1 be loops on the surface M representing the canonical generators 
of the fundamental group 7Ti(M). We require a,j to surround the points E 2 j_i, E 2 j 
(thereby changing sheets twice) and bj to surround E n , E 2 j-i counter-clock wise 
on the upper sheet, with pairwise intersection indices given by 

(2.4) di o cij — bi o bj — 0, ttjo bj — Sij, 1 < i,j < g. 

The corresponding canonical basis {Cj}j = i for the space of holomorphic differentials 
can be constructed by 

(") £ = t*)^. 

j = l K 2g+2 

where the constants c(.) are given by 

T^dn „ f E2k z^dz 



Cj(k) — C-^, Cjk — / 

J a 



pi/ 2 /„ R l/2 • 

The differentials fulfill 

(2-6) / Cfe = Sj,k, / (k=Tj,k, T jt k=T k ,j, l<j,k<g. 

Jaj Jbj 

Now pick g numbers (the Dirichlct eigenvalues) 
(2-7) &j) g j= i = fa, °i)U 

whose projections lie in the spectral gaps, that is, (ij G [E 2 j_i,E 2 j]. Associated 
with these numbers is the divisor T)^ which is one at the points (ij and zero else. 
Using this divisor we introduce 

(2.8) z(p,n)=A po (p)-a po (V^)-nA 00 _(oo + )-E po e C 9 , z(n) = z{^+,n), 
where S po is the vector of Riemann constants 

(2-9) %o, j = 1 ~ E l= lTj '\ Po = (E ,0), 

and A po (a po ) is Abel's map (for divisors). The hat indicates that we regard it as a 
(single- valued) map from M (the fundamental polygon associated with M) to C 9 . 
We recall that the function 8(z(p, n)) has precisely g zeros (ij{n) (with fij(0) = fij), 
where 6(z) is the Riemann theta function of M. 
Then our Jacobi operator H q is given by 

2 2 6(z(n + l))9(z(n-l)) 

a[U) ~ 8{z(n)f 

(2-10) b(n) = b + ±cM^4 J { - + r in) L ) ■ 

y ' K ' ^ ny 'dwj \6(w + z(n-l))J w=o 
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The constants a, b depend only on the Riemann surface and will be defined in the 
next section. 

It is well known that the spectrum of H q is purely absolutely continuous and 
consists of g + 1 bands 

(2-11) a(H q )=\J[E 2j ,E 2j+1 ]. 

3=0 

For further information and proofs we refer to [2E], Section 9. 

3. The Baker- Akhiezer function and the quasi-momentum map 

The Baker- Akhiezer function ip q {jp, n) = ip q {p, n, 0) is given by 
(3.1) 

\ 6{z{n a -l))0{z{^) 9(z(p,n)) / x 
n > «o) = W g(z(n _ imz(n)) 8(z( p ,n n )) 6XP K {n no) L ) ' 



po 



where Woo + ,oo_ is the normalized Abelian differential of the third kind with sim- 
ple poles at oo± and residues ±1, respectively. They are normalized such that 
ip q {p,n ,n ) = 1. 
The two branches 

n-l 

(3-2) ipq,±(z,n) = Y[ <t>q,±( z >j)> 

3=0 



where ([2Hj, (8.87)) 

(3.3) 0„ ± (z, n) = —- (z-b q (n)+Y: — ± 



2a g (n)\ '' j^z-fijin) Uj=i( z ~ ftW)/ ' 

#i W = n 7 7^ TTV » M = CT i W^j M » 



of the Baker-Akhiezer function are solutions of r g u = zh, z S C, where r q is the 
difference expression associated with H q . However, the Wronskian 

,1/2 

(3-4) </>,,+ (*)) 



n?=i(*-/^-) 

(/i 3 = /ij (0)) shows that they are linearly dependent at the band edge Ej, < j < 
2.9 + 1. 

The branch ip q ^ aj (z,n) has a first order pole at [ij if Hj is away from the band 
edges 

(3.5) lim (z - Hj)ipq, ffj {z, n) = tp q ^ (pj,n, l)-^r 

Z d q \ \) j 

(use ()3.3J) and ipq,±(z, n) = ip q ,± (z, n, l)4> qj ±(z, 0)) and both branches have a square 
root singularity if \ij coincides with a band edge Ei 



i' n M/ \j\Ei - E k \ 

(3-6) lim y/z- fj,jip q ,±(z,n) = ± - ip q +(Ei,n, 1) 

z ^ 2a 9 (0) n fe ^ V^i - Mfc 
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Lemma 3.1. The solutions of r q u — zu can be characterized as follows, 
(i) If R,2 g +2(z) 0, there exist two solutions satisfying 



( [ {z ' +) \ 

(3.7) ip q> ±{z,ri) = 8±(z,n)w(z) ±n , w(z) = exp ( / £>oo + ,oo_), 



with 9±{z,n) quasi-periodic. 

(ii) If R2g+2{z) = 0, z = Ei, there are two solutions satisfying 

(3.8) i/> q {E h n) = i/> q ,+ (E h n) =rl> q ,-(Ei,n), $ q (E h n) = il> q {E h n)@i(n) + n), 

where 0i(n) is quasi-periodic. 

Proof, (ii). We construct a second linearly independent solution at z = E = Ei 
using (see (1-50)) 



(3.9) s q (E,n) = lim a, (CI) 



V'g.+ fc.rc) - ip q -{z,n) 
ZTe"^"' Wty q ,-(z),ip q>+ (z)) ' 



where s q {z, n) denotes the fundamental solution of r q u = zu with initial conditions 
s q (z, 0) = 0, s g (z,l) = 1. W.l.o.g. we assume that Ei does not coincide with 
one of the Dirichlct eigenvalues fij (otherwise shift the base point). To derive an 
expression for ip q< ±(z) a,t z = E + e 2 we start with 

R 1 2 / g 2 + 2(z) = e(R + 0(e 2 )), R=-\{^E^E~. 

3+1 



Moreover, 

R_ 

and for p = (E + e 2 , ±) (see (|3.11ll below), 



W(ip q ,-(z),ip q ,+ (z)) = — g /jp ^ e(l + 0( £ 2 )) 



and 



^oo+,oo_ = / w±/3e + (9(e J ), /3= J - — g . 

Po Jpo R 

9 2E : >~ 1 

z(p,n) = z(E,n)±2e + 0(e 3 ), 7 = — S — ' 



<9t9 

U(P, »)) - 0(z(E, n)) ± + 0(e 3 ). 



Using this to evaluate the limit e — > shows 

s 

E - A: 



a g _ 

s q (E,n) = 2a q (0)\[—^^ q {E 1 n) = 1> q (E,n)0(n) +n) 



3=1 

where 



0(n) = t i ^)^ln^,») +ffi ), 
and finishes the proof. □ 
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Remark 3.2. (i). Since ip q (z,n) has a singularity if z = fj,j the solutions in 
Lemma \S.l\ are not well-defined for those z. However, you can either remove the 
singularities ofipq(z, n) or choose a different normalization point no ^ to see that 
solutions of the above type exist for every z. 

(ii). In the periodic case Floquet theory tells you that there are two possible cases at 
a band edge: Either two (linearly independent) periodic solutions or one periodic and 
one linearly growing solution. The above lemma shows that the first case happens 
if the corresponding gap is closed and the second if the gap is open. 

To understand the properties of ^±(2,71) we need to investigate the quasi- 
momentum map 

(3.10) uj(z)=expf / £Joo + ,oo_), p=(z,+). 

The differential w 00+j00 _ is given by 

/onx Il|=i(^-Aj) 

(3.11) Woo + ,oo_ = — — 572 dn, 

R 2g+2 

where the constants Xj have to be determined from the normalization 

(3.12) f = 2 / rc^/i*" X ^ dz = °> 

which shows Xj £ {Ezj-i, E2j). 

Since Xj G (i?2j-i, E%j) the integrand is a Herglotz function and admits the 
following representation (c.f. j^Hl, Appendix B) 

(3 - 13) < \ = / ~X~z MX} 

R 2g+2( Z ) J ~°° A Z 



with the probability measure 

7rii? 2g+2( A ) 



(3-14) dJkiX) = ll3 ,- 1 1/2 J^ XoiH^dX. 



Hence 

g(z,oo) = / tJoo + ,oo_ = / / x - / . djl(X)dC 

J pa J E • 

(3.15) = J™Jn(^f)d}l(\L 



A-C 



In particular, note that — Ke(g(z, 00)) is the Green's function of the upper sheet 
LT + with pole at 00+ and fj is the equilibrium measure of the spectrum (see |30| . 
Thm. III. 37). We will abbreviate g(z) = g(z, 00). 

The asymptotic expansion of exp(g(z)) is given by (|28|. (9.42)) 



(3.16) exp^ P ^ oo+:OO _^-^(l + ^ + 0( 1 



00, 
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where a is the capacity of the spectrum and 

1 2 ff +2 s 
(3-17) 6=^ E^-EV 

Theorem 3.3. TTie map g is a bijection from the upper (resp. lower) half plane 
C ± = j> S C | ±Im(z) > 0} to 

g 

(3.18) S ± = {z e C | ±Re(z) < 0,0 < Im(z) < tt}\ (J [ 5 (A,), g(E 2j+1 )\ 

3=1 

smc/i £/ia£ cr(H q ) — {z | Re(z) = 0}. 

Proof. By the Herglotz property of its integrand, the function oo) satisfies the 
conditions of [23], Theorem 1(b) in Chapter VI, which shows that it is one-to-one. 

To prove that g{z, oo) is surjective, it suffices to show that the boundary of C + 
is mapped to the boundary of S + . Note that g(X) is negative for A < Eq and purely 
imaginary for A S [Eq,Ei]. At E\, the real part starts to decrease from zero until 
it hits its minimum at Ai and increases again until it becomes at E2 (since all 
a-periods are zero), while the imaginary part remains constant. Proceeding like this 
we move along the boundary of S + as A moves along the real line. For A > E2 g +i, 
g(X) is again negative. □ 

Remark 3.4. In the special case where H q is periodic the quasi-momentum is given 
by w(z) = exp(iiV _1 arccos A(z)), where A(z) is the Floquet discriminant, and our 
result is due to pi] . 

Therefore the map 

9 

w:C ± -» W ± = {w e C I M < l,±Im(w;) > 0}\ \J[w(Xj),w(E 2j +i)] 

3=1 

(3.19) z 1— > exp(g(z)) 

is bijective. Denote W = W+ U W~ U (-1,1), W = W\{0}. If we identify 
corresponding points on the slits [w w (-©2^+1)] w e obtain a Riemann surface 
W which is isomorphic to the upper sheet II + . 

Remark 3.5. In |24) the largest band edge i?2g+i *s chosen for pq and w will map 
— > W T in this case. Moreover, in the periodic case the slits [w(Xj),w(E2j+i)] 
appear at equal angles where N is the period. 

Since z 1— > w(z) = exp(g(z)) is a bijection, we consider the functions ip„ ± as 
functions of the new parameter w whenever convenient. For notational simplicity 
we will write ip q ,±(w,n) for ipq,±(M w )i n ) an d similarly for other quantities. The 
functions ipq,± ( w : n ) are meromorphic in W and continuous up to the boundary with 
the only possible singularities at the images of the Dirichlet eigenvalues w(/j.j) and 
at 0. More precisely, denote by M± the sets of poles (and square root singularities 
if fj,j = Ei) of the Weyl m-functions m±{X), i.e. M + UM_ = {^j} 9 J=1 (see (j3.2jl and 
PS]. Section 2.1). Note that fij e M + n M_ if and only if fij = E t . Then 

(Bl) ip q ±(w,n) are holomorphic in W\({w(/j,j)}? =1 U {0}) and continuous on 
0W\{w(fj,j)}. 
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(B2) ifj q ,±(w, n) has a simple pole at w(fj,j) if fij € M±\{Ei}, no pole if fij ^ M±, 
and if fij = Ei, 

^,^=±^ + 0(1), 
w — Wl 

where C(n) is bounded and real. 
(B3) ip q ,±(w,n) =ij) q ^(w,n) for \w\ = 1. 
(B4) At w ~ the following asymptotics hold 

%M^n) = (-1)" ( ^W^ ) w ±»(l + OW). 



By Section 2.5 of [2H] the vector valued functions 
(3.20) U{X,n) 



1 f V><?,+ (A,n) 



4a g (0) 2 7rIm(m+(A)) V V-g,-(A,n) 

form an orthonormal basis for the Hilbcrt space L 2 (a(H q ),C 2 ,dX). The Weyl m- 
functions rh±(z) satisfy (see [21], ec t- (8-27)) 

(3.21) Im(m ± (A)) = t^^Y V A e o{H q ). 

r 0,+) 

thogonal basis on the unit circle 



Using our map w(z) = exp(f^ + Woo + ,oo_) we can transform this into an or- 



Lemma 3.6. Both functions ip q , + (w, n) and ijj q - (w, n) form orthonormal bases in 
the Hilbert space L 2 (5' 1 , ^dui), where 

,~ „„x , / \ ~f~T Mw) — fij dw 

3.22 duj(w) = H 7 —. 



Proof. Just use 
(3.23) 



2g+2\ 

□ 



Observe that is meromorphic on W with a simple pole at w = 0. In particular, 
there are no poles at w(Xj). 

Remark 3.7. In the periodic case we have 

N W_1 A-A 

(3.24) U p ..±W\\ 2 n := E l^iC^n)! 2 = N J] 1. 

n=l j=l A ™ 

4. Existence of Jost solutions 

After we have these preparations out of our way, we come to the study of short- 
range perturbations H of H q associated with sequences a, b satisfying a(n) — > 
a q (n) and b(n) — > as |n| — > oo. More precisely, we will make the following 

assumption throughout this paper. 
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Hypothesis H.4.1. Let H be a perturbation such that 

(4.1) H (l fl ( n ) " a »l + l 6 W ~ Wl) < oo. 

n£Z 

We first establish existence of Jost solutions, that is solutions of the perturbed 
operator which asymptotically look like the Baker- Akhiezer solutions. 

Theorem 4.2. Assume (H \4-l\ ). Then there exist solutions ip±(z,.), z € C, of 
Tip = zip satisfying 

(4.2) lim \w(z)* n (ip ± (z,n)-Tp q ,±(z,n))\=0, 

71 — >±OG 

where tp q ,±{z,.) are the Baker- Akhiezer functions. Moreover, ip±(z,.) are contin- 
uous (resp. holomorphic) with respect to z whenever tpq,±{ z i-) are an d inherit the 
properties (Bl) and (B2), where now tf)±(z,n) = 1 7 ±(n) + 0(1) and (54) has to 
be replaced by 

(4.3) xP±(z,n)=A±(0) 



nr= ; t a(m) ) (i+(B ± (o)±j2b(j-«))i+o(^ 

) 3 = 1 



where 



oc 



A+(n) = ft^y, B + (n)= £ (b q (m) - b(m)), 

j—7i ?n=n+l 

71 — 1 / .\ 71 — 1 

(4.4) A_(n) = [] ^> S -W= E (6 fl (m)-6H). 

j=— oo m=— oo 

Proof. The proof can be done as in the periodic case (see e.g., ^7], [23] or |28|. 
Section 7.5). The only problem is to show that the second solution at a band edge 
grows at most linearly. In the periodic case this follows from Floquet theory, here 
we just use Lemma f3. II □ 

From this result we obtain a complete characterization of the spectrum of H. 

Theorem 4.3. Assume (H \4-.l\j . Then we have o~ ess (H) = a(H q ), the point spec- 
trum of H is finite and confined to the spectral gaps of H q , that is, cr p (H) C 
M\o~(H q ). Furthermore, the essential spectrum of H is purely absolutely continu- 
ous. 

Proof. Again the proof can be done as in the periodic case (see e.g., [2H] or [2"£| . 
Section 7.5). □ 

5. The transformation operator 

We define the kernel of the transformation operator as the Fourier coefficients 
of the Jost solutions ip±(w,n) with respect to the orthonormal system given in 
LemmaED {ip q ,±(w,n)} neZ , 

(5.1) K±{n,m) := / ip±(w, n)ip q . T (w, m)dui(w). 

27TlJ| TO | =1 

By the Cauchy theorem, this integral equals the residue at w = 0, 

(5.2) K±(n, to) = Reso— tp±(w, n)ip q T (w,m). 

w 
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In particular, since ip±(w , n)ip q . T (w , m) = 0(w ± ^ n ^ m ^), we conclude 

(5.3) K±(n,m) = Q, ±(m-n)<0. 

Lemma 5.1. Assume H \4- 1\ The Jost solutions tp±(w,n) can be represented as 

±00 

(5.4) i(;±(w,n) = ^ K±(n,m)ip qy ±(w,m), \w\ = 1, 

m—n 

where the kernels K±(n, .) satisfy K±(n,m) — for ±ra < zbn and 

±00 

(5.5) \K ± (n,m)\<C £ (|a(j) - a,tf)l + - ^COl), ±m > ±n. 

j=[^]±i 

TTie constant C depends only on H q and the value of the sura in 

Proof. We prove the estimate for K + (n,m) and omit "+" and "z" whenever pos- 
sible. Define = ip(n)K(n, n) -1 , then fulfills 

00 

(5.6) (p(n) — ip q (n) + J(n,m)ip(m), 

m— n+l 

where 

/r - x, x -/ ^ s q (z,n,m- 1) ~, s g (z,n,m) 

(5.7) J(z,n, m) = aim — 11 ; — V b{m) — — ; — - — 

a q (m — 1) a q {m) 

with the abbreviation 

(5.8) a(m) — — - — - — a q (m), b(m) = b(m) — b q (m). 

a q (m) 

On the other hand, <p(n) is given by 

f( n ) = y] K(n,m)t/) q (m), n{n,m) = ^J 1 '™, , 

K(n,n) 

m=n v 7 

therefore 
(5.9) 



CO oo 



/t(n, m)ip g (m) = J{n,m)^p q (m) + J(n,m)K(m,l)tp q (l). 

m—n m— n+l m— n+l /— m+1 

Multiplying both sides of (|5.9|) by (fc) and integrating over the unit circle yields 

oo oo oo 

(5.10) K,(n, k) = r(n, m, to, fc) + T(n,m,l,k)K,(m,l), 

m=n-\-l Tci— n+l i=n+l 

where 

(5.11) r(n,m,/,fe) = / J(w,n,m)ip q + (w,l)ip q -(w, k)dui(w). 

2tti J| w |=i 

Using EU, (1.50), 

, g s q (n,m) = ip q ,+ (m)ip q -(n) - V> g , + (n)-0 g _ (to) 

flt(m) W(V>?,+ ,^g,-) 

we obtain 

(5.13) r(n, to, Z, fc) = b(m)T q (n, to, Z, fc) + a(m)r g (n, m — 1,1, k) 
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with 



T q (n, to, I, fc) = Ta(m,n,l,k) — Ta(n,m,l,k), 
_j_ f ^ g ,+ (w,n)ip q -(w,i 

1 f , 4'q,+ ( Z , n )^q.~( Z ,m)lp qt+ (z,l)lp q ^(z,k)Yl(z - flj) 



p / t ; \ 1 f i>q,+ {w,n)^q-(w,rn)'ip q<+ (w,l)'ip q -(w,k) 

ln(n,m,l,k) — — : / — — 77- ; — : — -, ; — rr— dioiw) 



2m Wty q ,+ (z) t %,-{z)) R\' 2 +2 {z) 



dz 



r r U s 1 / V'g.+ fo n)ljj q -(z, m^+jz, l)4>q-{z, fc) 

1 J ~ 2^i7 7 W(^, + («),V g> -(«)) 2 

Here 7 is a path on the upper sheet encircling the spectrum. The integrand of To is 
meromorphic on the Riemann surface M with poles of order one at Ej and poles of 
order O( ;s ± 0»- m +*- fc )- 2 ) near oo± (there are no poles at the Dirichlet eigenvalues 
fij ). We apply the residue theorem twice, first on the side of 7 including oo + , then 
on the other side including the spectrum (and thus 00 _) 



r (n, m, Z, fc) = —Res, 



V> g , + (n)V>,,- (rn)if)g }+ (l)ip q ,- (fc) 



2g+l 



(5.15) = (Rc Soo _ + Res Bj 

3=0 

The order of the poles at oo± implies 

2 S +1 



(")^,-( m )^,+ (0 V'g - ( fc ) 

W^+,V 9 ,-) 2 



T (n,m,l,k) = \ ,e J w^T^? n m + / fc<0 

n-m + Z- fc>0, 

which shows that To(n,m,l,k) is real and bounded since ij) q ^{E,.) = ip q ^(E,.) 
are (if [i 3 = Ei, use (B2)). Together with (|5.14(l this yields 



(5.16) ro(n, m, Z, fc) = — ro(m, n, fc, Z) = — ro(m, n, fc, Z) = — To(n, m, fc, Z). 
Moreover, 

Y q (n,m,l,k) = 0, Z— fc > |m — n|, 

(5.17) r g (n, m, Z,fc) = — r 9 (m, n, fc, Z) = T q (n,m,k,l), 

which then implies 
(5.18) 

r>, m , z, fc) = ( si ^ - ^)g W- +(w ^^r'- (,t) I' - *l < l™ - n l 

[ |Z - fc| > |to - n| 

and r(n, to, Z, fc) = for |Z — fc| > m — n if m > n. Note that the residue at Ej is 
given by 

(5.19) 2 ffH ( f J '~p| ^ ' ^ ' ' ' fc) ' 
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Now we obtain for ^(n, k) 

OO 

n(n 7 k) = r(n, m, m, /c) + T(n 1 rn 1 l 1 k)n{rn 1 l) 



CO CO 



m— n+1 



m— n+1 /— m+1 

oo m+fc— n— 1 



(5.20) 



r(n, m, m, fe) + T(n,m,l,k)K,(m,l), 



m— n+1 l—n+k—m+1 



since T(n,m,m, k) ^ only if |m — fc| < m — n implying m > ^j^. In the third 
sum of 15.20fl we need that |m + <5 — fc| < m — n for <5 > 1 which yields 5 < k — n 
and <5 > n + fc — 2m. Two remarks might be in order: m+k—n—l > n + k — m + l 
since m — n > n — to + 2, and the starting point / = n + fc — m + 1 of the third 
sum actually has a lower limit, namely m < ^y 1 , since we require I > to + 1 for 
«(m,Z) ^0,1. Note that 

oo 

r(n, m, m, k) 



m =[2jS]+l 

m+A;— n — 1 



< D J] |6M+a(m)|=:$(a^), 

ro= [2±*]+l 



X] |r(n,m,Z,fc)| 



< D(m-n - l)|6(m) + a(m)\ =: c(m) € ^(Z), 

Z— n+fc— m+l 

where I? is the estimate provided by Ij5.18|l , (|5.19(1 . We set up the following iteration 
procedure 



Ko(n, k) 



^2 T(n,m,m,k), 

m= [B±*]+l 

oo m+/c— n— 1 

(5.21) Kj(n,k) — r(n, m, i, fc)Kj_i(m, Z). 

m=n+l /— n+fe— m+1 

Then using induction one has 



(5.22) 



i«i(M)i<«m 



n+ fe,(Em=«+l 2 ( m ))' 



and hence the iteration converges and implies the estimate 



(5.23) 



|/t(n, A;) | 



^Kj(n,fc) 



<9(^)cxp ^ c(m) 



i m— n+1 



□ 



Associated with K±(n, m) is the operator 



±oo 



(5.24) 



(£±/)M = J2 K ± (n,m)f(m), f G *£(Z,C)> 



which acts as a transformation operator for the pair r, r q . 

Theorem 5.2. Let T q and r be the quasi-periodic and perturbed Jacobi difference 
expression, respectively. Then 



(5.25) 



t/C±/ = /C ± V, /e*£(Z,C). 
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Proof. It suffices to show that HK± = K±H q . 



HK±(n,m) = / HiJ;±(w,n)ip q . T (w,m)duj(w) 

\w\ — l 



2ni 
1 
2tti 



A(w)ip± (w, n) , gj =p(iy, m)duj(w) 



w\ — 1 



(5.26) = / ijj±(w,n)H q ip qtT (w,m)duj(w). 



2ni 



\w\=l 



□ 



Lemma 5.3. For n E Z we have 

a(n) K + (n + l,n+l) K_(n,n) 



(5.27) 



a q {n) K + (n,n) K-(n + l,n+ 1) 

if + (n,n + l) , , 2f+(n-l,n) 
b(n)-b q (n) = a q (n) — — — r a q (n-l)- 



K + (n,n) qy ' K+(n-l,n-l) 



K-(n,n) 9W if_(n + l,n+l)' 

Proof. Consider the equation of the transformation operator HJC± = K,±H q , which 
is equivalent to (c.f. I|5.26[) ) 

a(n — l)K±(n — 1, m) + b{n)K±(n, m) + a(n)K±(n + 1, m) = 

= a q (m — l)K±(n, m — 1) + b q (m)K±{n, m) + a q (m)K±(n, m + 1). 
Evaluating at m = n we obtain the first equation and at m = n =p 1 the second. □ 

In particular, observe 

(5.28) K ± (n,n) = A±(n), K±(n,n± 1) = -^%-B±(n). 

a q {n- J 



6. The scattering matrix 

Let -ffq be a given quasi-periodic Jacobi operator and H a perturbation of H q 
satisfying Hypothesis H I4.1I To set up scattering theory for the pair (H,H q ) we 
proceed as usual. 

The Wronskian of our Jost functions can be evaluated asn-> ±oo and is given 

by 
(6-1) 



W(<MA),V±(A)) = W q (% t± {X),% tT (X)) = T n9 ^ v A e a ( H i 



Hence V'±(A), ^±(A) are linearly independent for A in the interior of o~(H q ) and we 
consider the scattering relations 



(6.2) V±(A,n) =a(A)^ T (A,n)+/3 T (A)^ T (A,n), A € cr(H q ), 
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where 

(6 . 3 , o(a, = «44§ = r%^i) 

W / (V'±(A),V±(A)) fl&+ 2 (A) 

While a(A) is only defined for A S a(H g ), 1)6.3(1 may be used as a definition for 
A e C\{£j}. Therefore can be continued as a holomorphic function on W 

and it is continuous up to the boundary except possibly at the band edges. 

Remark 6.1. Note that a(X) does not depend on the normalization of ip±(X) at the 
base point n$ — whereas (3 ± = 0± t o does. Using ip±(z,n, no) = ipq,±(z,no)~ lr ip±(z 1 n) 
and 

W(ty+(A),V-(A)) = II X ~^ j{no) W((i> + (\,.,no),i>-(\,.,no)) 

we see 

(6-4) P ± AX) = ^ n °W n (X). 

A direct calculation shows 



(6.5) a(w) = a(w), (3±(w) — (3±{w) = —(3^{w) 
and the Pliicker identity (c.f. [2H|> (2.169)) implies 

(6.6) \a(w)\ 2 = l + \(3 ± (w)\ 2 , H = l. 
We will denote the eigenvalues of H by 

(6.7) a p (H) = { Pj } q j=v 

Our next aim is to study the behavior of a(A) at the eigenvalues pj, therefore we 
modify the Jost solutions ip±(X,n) according to their poles at fij and define the 
following cigenfunctions tp±(X, .) 

(6.8) i£+(A,.)= [] (A-w)lMV), ^-(A,-)= II (A-w)^-(A,.)- 

Define t/> g ,±(A, .) accordingly. Moreover, 4>±{pj,n) = c^ip^(pj,n) with c^cj = 1. 
The norming constants 7±.j are defined by 

(6.9) J_ := J2\4>±(Pi,m)\ 2 ' 
To compute the derivative of a(X) at pj, note that 

(6.10) a(A) = 1/2 . 

(A) 

By virtue of Lemma 2.4, 

(6.11) A^(VL(A),^ + (A)) = -^_( ft ,fc)^ + ( ft ,fc) 
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Therefore 

(6.12) *a{X) - U " lV - ' ± 



From (|6~l"2"|) we obtain a connection between the left and right norming constants 

( 6 - 13 ) 7+j7-j = - f n vJb 7 — \- 

As a last preparation, we study the behavior of a(w) as w — > 0. By l|4.3|) . 

(6.14) W(^_(V), ^+(to)) - ^(OjA+^ouT 1 + O(w) 

and 

therefore a _1 (w) is bounded at with 

(6.16) «(0) = g ^ 

j'=— oo ^ ' 

We now define the scattering matrix 

(6.17) 5(t») = 



T{w) R-(w) 
R+(w) T(w) 

where T(w) :— a^ 1 (w) and R±(w) :— a~ 1 (w)[3±(w) are called transmission and 
reflection coefficients. Equations (|6.5(l and 1)6.6(1 imply 

Lemma 6.2. The scattering matrix S(w) is unitary. The coefficients T(w), R±(w) 
are bounded for \w\ — 1, continuous for \w\ — 1 except at possibly wi = w(Ei), fulfill 

(6.18) \T{w)\ 2 + \R ± {w)\ 2 = 1, H=l, 

(6.19) T{w)R+{w) +T(w)R-(w) = 0, \w\ = 1 



and T(w) = T(uJ), i?±H = i?±(w) /or |io| = 1 

Moreover, R\^ + 2{ w )T{w) _1 is continuous (in particular T(w) can only vanish 
at wi ) and 

lim R\ / g 2 +2 iw) R± ^) 1 = °. ^ ^ ">(^) 

^ ' V r>l/2 / x.R ± (u>)-l n / \ " 

J™, ^23+2 ( w ) T(J) = °« ">» = W W> 

The transmission coefficient T(w) has a meromorphic continuation to W with 
simple poles at w{pj), 

(6-21) (Res Pj T(A)) 2 = 7+ , i7 _ >j ii 2ff+2 ( (0j )- 

In addition, T{z) £ M as z 6 R\cr(iJ g ) and 

^ r (0) = 1 r = fl 44 

K + ti.n «_ ti.n ■ L1 - a„(m 

\ ' / m— — oo y v ' 

where K±(n,n) are the coefficients of the transformation operators. 
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Proof. To show (I6.20|) we use the definition H6.3(l . 

rtiU ^TiX) 1 =]I( A -^)(^(V'-(A)^ + (A))T^(^(A),^(A))). 

There are two cases to distinguish: If [ij ^ Ei then ip± are continuous and real 
at A = Ei and the two Wronskians cancel. Otherwise, if /ij = Ei they are purely 
imaginary (by property (B2) of the Jost functions) and the two terms are equal in 
the limit and add up. □ 

The sets 

(6.23) S ± {H)={R ± (w),\w\ = l; ( Pj ,7 ±tj ), 1 < j < q} 

are called left/right scattering data for H. 

First we want to show that the transmission coefficient can be reconstructed 
from either left or right scattering data. 

Let g(w,wo) be the Green function associated with W and let 

(6.24) Li(w,wo)dw = ^-(w,re ie ) e w d9, w = e w , 

or r=l~ 

be the corresponding harmonic measure on the boundary (see, e.g., Since 
Wo is simply connected, we can choose a function h(w,v) such that g(w,wo) — 
g(w,wo) + ih(w,wo) is analytic in Wo- Clearly g is only well-defined up to an 
imaginary constant and it will not be analytic on W\{0} in general. Similarly we 
can find a corresponding v(w, wo) and set fi(w,wo) = /z(u>,wo) + iv{w, wq). 

Theorem 6.3. Either one of the sets S±(H) determines the other and T(w) via 
the Poisson- Jensen type formula 
(6.25) 



T(w) - exp (^gfawfa))^ exp ^ln(l - \R ± (i 



w )\ )fi(w,wo)dwo 

where the constant of g has to be chosen such that T(0) > 0, and 
R-(w) T{w) (Res Pj .r(A)) 2 
R+(w) T(w) ' 7+j7 - J nf^fo - E,) ' 

Proof. It suffices to prove the formula for T(w), since evaluating the residua pro- 
vides 7±.j, together with {A;}, {£"/}• The formula for T(w) holds by [35], Theorem 1 
at least when taking absolute values. Since both sides are analytic, and have equal 
absolute values, they can only differ by a constant of absolute value one. But both 
sides are positive at w = and hence this constant is one. □ 

Note that neither the Blaschke factors nor the outer function in l|6.25|l are single 
valued on W in general. In particular, the eigenvalues cannot be chosen arbitrarily, 
which was first observed in |21| . 

7. The Gel'fand-Levitan-Marchenko equations 

In this section we want to derive a procedure which allows the reconstruction 
of the Jacobi operator H with asymptotically quasi-periodic coefficients from its 
scattering data S±(H). This will be achieved by deriving an equation for K±{n, m) 
which is generally known as Gel'fand-Levitan-Marchenko equation. 
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Since K±(n,m) are essentially the Fourier coefficients of the Jost solutions 
4>±{w, n) we compute the Fourier coefficients of the scattering relations ()6.2JI . There- 
fore we multiply 



(7.1) T(w) , ipzp(w, n) = R±(w)ip±(w, n) + ip±(w, n) 

by (2iri)~ 1 ift qt ±(w 1 m)duj, where ±m > ±n, and integrate around the unit circle. 
First we evaluate the right hand side of (|7.1|) using (|5.1|l 



(7.2) — / ip + (w,n)?p q ^ + (w,m)du(w) = K + (n,m) 



1 

27ri 



w\ = l 



oo 

R + (w)^J + (w,n)ip q .+ {w,m)duj(w) — K + (n,l)F + (l,m), 



l—n 



where 



(7.3) F + (l,m) = — I R + (w)i/> q:+ (wJ)^ qt+ (w,m)dw(w) 

\w\ = l 



2ni 



Note that F + (l,m) = F+(m,l) is real. 

To evaluate the left hand side of (|7.1|) we use the residue theorem. The only 
poles are at the eigenvalues and at if n = m, hence 



— — : / T(w)ip^(w,n)ip g + (w,m)duj(w) 

2ln .l\ w \ = l 



S(n,m) 



K + (n,n) «V ^ a(A) 

Here S(n, m) is one for m = n and zero else. By (|6.12|) the residua at the eigenvalues 
are given by 

(7.4) Res, «f X > m) ) = - 7+ ^ + ( ft ,n)^ + (p,,m). 

Collecting all terms yields 
(7.5) 

±°o ^, s q 

K±(n,m) + y2K±(n,l)F ± (l,m) = c -V7iii(p„n)ii(ft,H 

^— ' K±(n,n) ' 

and we have thus proved the following result. 

Theorem 7.1. The kernel K±(n,m) of the transformation operator satisfies the 
Gel'fand-Levitan-Marchenko equation 



(7.6) K ± {n, m) + y K±(n, l)F±(l, m) = 5"/ ™\ , ±m>±n, 



wh 



ere 



(7.7) F ± (l,m) = F ± {l,m) + ^- 1± ^ qt± { Pv l)^ q ^{ Pj ,m) 

3=1 
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Defining the Gel'fand-Levitan-Marchenko operator 

oc 

(7.8) ^±/(j) = £i^(n±Z,n±j)/(0, /e^ 2 (N ,C), 

1=0 

yields that the Gel'fand-Levitan-Marchenko equation is equal to 

(7.9) (l+J*)K ± {n,n±.) = (^(n.n))" 1 ^. 

Our next aim is to study the Gel'fand-Levitan-Marchenko operator in more 
detail. The structure of the Gel'fand-Levitan-Marchenko equation suggests that 
the estimate 1(5. 5() for K±(n 7 m) should imply a similar estimate for F (n,m). 

Lemma 7.2. 

±oo 

(7.10) \F±(n,m)\ <C £ (|a(j) - a q (j)\ + \b(j) - b q (j)\), 

where the constant C is of the same nature as in J^5.5}) . 
Proof. We abbreviate the estimate (|5.5ll for K + (n,m) by 

(7.11) \K + (n,m)\<CC + (n + m), 
where 

oo 

C + (n + m)= C W> c(j) = \a(j)-a q (j)\ + \b(j)-b (! U)\. 

Note that C+(ra + 1) < C+(n). Moreover, C+(n) £ ^+(Z) since the summation by 

parts formula (e.g. [28: , (1.18)) 

(7.12) 

N N 

<?(m)(/(m+l)-/(m)) = g(N)f(N+l)-g(n-l)f(n)+J2(9(m-l)-g(m))f(m) 

m—n m—n 

implies for g{m) — m. f(m) — C+(m) that 

OO CO 

(7.13) £mc(m) = (n-l)C + (n)+£c + (m), 

m—n m— n 

where we used linin^+oo n C+ (n + 1) < linin^oo J2m=n m c ( m ) — 0. Solving the 
GLM-equation 1)7.6(1 for _F + (n, m), m > n, we obtain 

\F+(n,m)\ < - ] l\K+(n,m)\+ V \K+(n, l)F+(l, m)\ ) 
< d(n) lc+(n + m) + C+(n + I) \F+(l, m) | J , 
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where Ci(n) = C \K+(n, n)^ 1 — ► C for n — ► oo (see l|5.28|l ). For n large enough, i.e. 
Ci(n)C+(2n) < 1, we apply the discrete Gronwall-type inequality [2Hjj Lemma 10.8, 

F+(n,m < d(n) C + (n + m) + ^ 1W + v / + v 

(7,4) < Ot(n)C + (» + m) fl+ £ , ^ffi±^ ) , 

V iMi nU+i(l - Ci(n)C + (n + *)),/ 

which finishes the proof. □ 
Furthermore, 

Lemma 7.3. Let F^{n,m) be solutions of the Gel'fand-Levitan-Marchenko equa- 
tion. Then 

(7.15) M \ F± ( n > n ) ~ F ± (n± l,n± 1)| < oo, 

n— no 

(7.16) ^ |n| |<Zg(n)F ± (n,rc+ 1) - a,(n - LjF^n- l,n)| < 00. 



Proof. We first prove (|7.16|) for F + . Lemma 15 . 31 implies 

(7.17) b(n) - b q (n) = a q {n)n + s{n) - a q (n - l)/c+,i(n - 1), 

where 

K + (n,n + j) 



(7.18) k +j (n) := K+(n, n + j) 



K+{n,n) 



Abbreviate iy(n) := F + (n + j,n). With this notation, the GLM-equation (|7.t)|) 
reads 

(7.19) K+ ,(n) + F+(n) + £ ,(n + /) - / ( /'°\ 2 , J > 0. 

Insert the GLM-equation for F + (n, n+1), F + (n~l, n) (recall F + (n, m) = F + (m,n)) 

a q {n)F+{n) - a q (n - 1)1? (n - 1) 

= -a q (n)K + , 1 (n) + a q (n — l)/t +) i(n - 1) 

00 

(7.20) - («>)k +j (n)F+ > + 1) - a q (n- 1)« +J (n - lJ^W) . 

3=1 
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Since — a g (n)tt+,i(ra) + a q (n — 1)k+,i (n — 1) = b q (n) — b(n) the only interesting part 
is the sum. For N, J < oo, 

N J 

^2 { a g(n)K+,j(n)Ff_ 1 (n + 1) - a q (n - l)«+j(n - ^n)) 

n— no J — 1 

,7 JV 

= X X ra(a«(ra)«+,j(")- F jti(n + 1) - - l)t+,j(n - l^t^n)) 

j — 1 n— no 
,7 

= £ (^W«+,iW^,ti(iV + 1) - (n - l)a 9 (no - l)K+,i(no - ljf+^W 
3=1 

7Y 

(7.21) + J2 (-^(n-lKitn-l)^^)), 

n— no 

where we used the summation by parts. Estimates <|7.11l) . (|7.14f) imply for the first 
summand 

J J 
^Na^K+jWFf^N + l^ < J2\ N M N ) dc +( 2N + 3) c +( 2N + j + 1 ) 

3=1 3 = 1 

< \N\a q (N)CC+(2N+l), 
which holds uniformly in J, and (compare (|7.13|) ) 

(7.22) lim NaJN)CC+(2N + 1) = 0. 

TV— >oo 

Moreover, 

J N 

j=l n=n 

J JV 

- jv 1 ]™^^ 5Z \ a i( n - 1 ) K +,j( n - tfFj'-iin) 

j=l n— no 

oo oo 

< a q {n - l)CC+(2n + j)C+{2n + j + 1) < oo. 

j — 1 n— no 

Therefore |n||a g (ra).F + (n, n + 1) - a q (n - l)F + (n - l,n)| G as desired. To 

apply Lemma |5.3l for y7~ use the symmetry property F~(n,m) = F~(m,n). For 
B7.15JI . inserting the GLM-equation yields 

F + (n, n) - F + {n + l,n+ 1) = K~ 2 (n, n) - K~ 2 {n + 1, n + 1) 

OO 

+ X (« +J (n + l)F/(n + 1) - n +>j {n)F+{r 

3 = 1 

By dSHHJ), 

i T , 9/ n „ 9/ , ,i \a(n) + a„(n)\ -rr o,(j) 2 . , . ... 

\Kf(n,n)-Kf(n+l,n+l)\ < > J K J] _^-| (n) - o,(n)| 

a\n) j=n+1 a q {j) 

(7.23) < C\a(n) - a q (n)\, 

and the same considerations as above imply l|7.15|l . □ 
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Remark 7.4. The Gel'fand-Levitan-Marchenko equation is symmetric in K±(n, m) 
and F^(n,m), therefore we can invert the analysis done in Lemma \7.'J\ and obtain 
estimates for K±(n, m) starting with an analogue of estimate J 7. 1 Uj) for F^(n, m) 
and the estimates \7.15\) , j7.16\ ) (c.f Lemma \8.1\) . 

Theorem 7.5. For n G Z, the Gel'fand-Levitan-Marchenko operator : £ 2 — > £ 2 
is Hilbert- Schmidt. Moreover, 1 + is positive and hence invertible. 

In particular, the Gel'fand-Levitan-Marchenko equation \7.9jj has a unique solu- 
tion and S+(H) or S-(H) uniquely determine H. 

Proof. That F^ is Hilbert-Schmidt is a straight-forward consequence of our esti- 
mate Lemma 17.21 

Let / € ^ 2 (No) be real (which is no restriction since F + (n, I) is real and the real 
and imaginary part of H7.24|l could be treated separately) and abbreviate f n (w) = 

OO oo oo 

E f(i)^fti) = E /0') E F+ ( n +i> n + o/co 

3=0 3=0 1=0 

If °° 

= —. I R+{w) E f(j)^q,+ (w,n + j)ip q , + (w,n + l)f(l) du(w) 
2m /hl=i 

<? OO 

+ E E fUh+,k$q,+ (Pk,n + j)-4> q>+ (p k ,n + l)f(l) 

fc=l j",Z=0 

/ R+(w)f n (w)f n (w)duj(w) + ^2l+,k\fn(pk)\ 
'M=i fe=1 

(7.24) = 7? + H|/„H| 2 ^H+E^I/»(^)| 2 ' 
27rl /|«|=i ^ 



1 

2vr7 



2 



where R+(w) = R + {w)f n (w)(f n (w)) with = \R+{w)\ and /„(«;) = 

X^=o f(J)' l l ; q,+ ( w i n + i)- The integral over the imaginary part vanishes since 
-R+(u>) = R+(w) and we replace the real part by 

Re(i? + M) = + ,R + H| 2 - 1 - |i?+H| 2 ) = + i?+M| 2 + |TH| 2 ) - 1, 
(recall |i? + H| 2 + |T( W )| 2 = 1). This yields using £ |/(j)| 2 = aS¥ Jj«[=x \UH\ 2 duj 

oo g 

E7w(i+^)/o-) = £7+,*i/»(p*)i 2 

3=0 fc=l 



(7.25) 



~ J w =i (|i + ^+MI 2 + |TH| 2 )|/„H| 2 ^(r 



which establishes l + .7 r + > 0. According to Lemma lrT2*l |T(uj)| 2 > a.e., therefore 
— 1 is not an eigenvalue and 1 + F^ > e n for some e n > 0. □ 

To finish the direct scattering step for the Jacobi operator H with asymptoti- 
cally quasi-periodic coefficients we summarize the properties of the scattering data 
S±(H). 
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Hypothesis H. 7.6. The scattering data 

(7.26) S ± (H) = {R±(w),\w\=l;(p j , 1± , j ),l<j<q} 

satisfy the following conditions. 

(i). The reflection coefficients R±(w) are continuous except possibly at Wi 
w(Ei) and fulfill 



(7.27) R±(w)=R±(w). 
Moreover, \R±(w)\ < 1 for w ^ wi and 

2g+i 

(7.28) 1- \R±(w)\ 2 > C J] \w-Wi\ 2 . 

1=0 

The Fourier coefficients 

(7.29) F ± {l,m) = ^-.( R±(w)ij q .±{w,l)i/j q , ± {w,m)du J (w) 

2vri J\ w \ =1 

satisfy 



^(n.m)!^ «(?)» ?0')>0, \j\q(j)ef(Z) 7 

j=n-\-m 

^ |n||F ± (n,n)-F ± (n±l,n±l)| < oo, 

n— no 

|n||a 9 (n)i^ ± (n, n + 1) — a q {n — l)_F ± (n — l,n)| < oo. 

n— no 

('nj. TTie values pj £ M\a(H q ), 1 < j < q, are distinct and the norming constants 
7±,j> 1 < J < <Z,. a?"e positive. 

(Hi). T(w) defined via equation \b\ 25}) extends to a single valued function on W 
(i.e., it has equal values on the corresponding slits). 

(iv). Transmission and reflection coefficients satisfy and satisfies 



(7.30) 



lim (w - Wi) fl T ( rt\ +1 = °' w i ^ w (t J 'j) 
lim (w-wi) R y^Jr 1 = wi=w(fij) 

and the consistency conditions 



2 



R-(w) _ T(w) (Res Pj T(A)) 
R + (w) T(w) ' 7+J 7 - J " n ^\ pj El ) ■ 

Remark 7.7. Note that \7.28\j implies that ln(l — |i?±(w)| 2 ) is integrable and 
ensures that \6. 25}) is well-defined, at least as a multi valued function. Condition 
(Hi), which is void in the constant background case, shows that the the reflection 
coefficient and eigenvalues cannot be chosen independent of each other. 
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8. Inverse scattering theory 

In this section we want to invert the process of scattering theory, that is, we 
want to reconstruct the operator H from a given set S± and a given quasi-periodic 
Jacobi operator H q . 

If S± (satisfying H I7.6I (i)— (ii)) and H q are known, we can construct F ± (i,m) 
via formula 1|7.7|) and thus derive the Gel'fand-Levitan-Marchenko equation, which 
has a unique solution by Theorem l7.5l This solution 

K±(n,n) = (6 ,(1+F±)- X 6 ) X/2 

(8.1) K ± (n,n±j) = } r (6j, (1 + 

K±(n, n) 

is the kernel of the transformation operator. Since 1 + is positive, K±(n,n) is 
positive and we can set in accordance with Lemma 15. HI 

, , . „ K+(n + 1, n + 1) 

(8.2) o+(n) = a q (n)- ^ 



a_(n) = a q (n) 



K+(n,n) 
K-(n, n) 



K-(n+l,n+l)' 



, , x , , , . , ,K + (n,n + l) K+(n-l,n) 
b + (n) = b q (n) + a q (n) — — — a q {n-l)- 



K + (n,n) qy ' K + (n-l,n-l)' 

, / \ , , s , s K-(n,n — V) , . K-(n+l,n) 
6_(n) = b q (n)+a q {n-l)— ^ -L - a q (n)- 



K-(n,ri) y ' K-(n + l,n+ 1) 
Let H+, H- be the associated Jacobi operators. 

Lemma 8.1. Suppose a given set S± satisfies H \ 7. b\ (i)-(ii). Then the sequences 
defined in hH.ty) satisfy n\a±(ri) — a q (n)\, n\b±(n) — b q (n)\ G ^±(N). 

Moreover, ip±(X,n) = ^2^^ n K±(n,m)ip q: ±(X,m) , where K±(n,m) is the solu- 
tion of the Gel'fand-Levitan-Marchenko equation, satisfies T±ip± = Xip±. 

Proof. We only prove the statements for the "+" case. Define F + (n,m) by (c.f. 

O) 

q 

F + (l,m) = F + (l,m) +^i +J 4> q , + {p j ,l)'4)g, + (p j ,m). 

3=1 

Hypothesis H I7.(jl fi) implies 

oo 

(8.3) \F+(n,m)\<C £ q(j) =: C+(n + m), 

j—n-\-m 

oo 

(8.4) \n\\F + (n,n)-F + (n+l,n + l)\ < oo, 

n—riQ 
oo 

(8.5) Y H\a q (n)F + (n,n+l)-a q (n-l)F + ( n — 1, n)\ < oo, 

n— no 

since ipq,+ (Pj, n ) decay exponentially as n — > oo and J^j l+.j^Pq-.+ iPji O^g.+ CPi) •) 
form a telescopic sum. Note that C + (n + 1) < C + (n). 
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Set K+(n,m) := K+(n,m)K+(n,n) . Then as in the proof of Lemma [7.21 we 
obtain 

(8.6) \K+(n,m)\ < C+(n + m)(l + 0(1)). 

Now we have all estimates at our disposal to prove n\b+(n) — b q (n)\ 6 ^ 1 (N). By 
definition (c.f. JH2J), 

(8.7) b+{n) — b q {n) = a q (n)K+(n, n + 1) — a q (n — \)n + (n — 1, n). 

We insert the GLM-equation for K+(n, n + 1), n + (n — 1, n) and use estimate 1)8. 5p . 
the summation by parts formula, and estimates (|8.3() . (|8.6() in the same way as in 
Lemma IV.3I Similarly using (|8.4|l we see 

1 1 



n— no 



K\{n,n) ^(n + l,n + l) 



< oo. 



Equation (|8.2|l yields 
1 



K\{n,n) K\{n+l,n+l) 



a q {n) 2 V j .j n ^ 1 a 9 (j) 2 

The product converges and therefore |n||a+(n) 2 — a q (n) 2 \ S £ X (N). 
Next we consider ip + (\,n). Abbreviate 

(8.9) (AK + )(n 7 to) — a q (n — l)n + (n — 1, to) + a+(n)a' q 1 (ri)K + (n + 1, to) 

—a q (m — 1)k_|_(ti, m — 1) — a q (m)n + (n, m + 1) + (6+(n) — & g (m))/c+(n, m). 

Ai^ + = is equivalent to the operator equality H + K + = K + H q , which in turn 
implies that ^>+(A, n) satisfies H + ^p + — Xip+ 

(8.10) H+iP+ = H+K + %. + = K+H q ^+ = K+\% t + = XK+^+ = A^+. 
To show that AK + = we insert the GLM-equation into l|8.9|l and obtain 

oo 

(8.11) (AK+)(n,m)+ (AK + )(n,l)F+(l , to) = 0, m > n + 1. 

2=n+l 

In the calculations we used 

a q (n - l)F+(ra - 1, to) + b q (n)F + (n, m) + a q {n)F + (n + 1, to) = 

a q (m — l)F + (n, m — 1) + b q (m)F + (n, to) + a q {m)F + (n, m + 1) 

which follows from 17. 7|) . By Theorem 17.51 equation 18.11|l has only the trivial 
solution AK + = and hence the proof is complete. □ 

Now we can prove the main result of this section. 

Theorem 8.2. Hypothesis H \71] is necessary and sufficient for a sets S± to be the 
left/right scattering data of a unique Jacobi operator H associated with sequences 
a, b satisfying H \4-l\ 

Proof. Necessity has been established in the previous section. By Lemma f8. II we 
know existence of sequences a± , b± and corresponding solutions {w, n) associated 
with S+ (or S-). Hence it remains to establish a + {n) — a_(n) and b+(n) = b-(n). 
Consider the following part of the GLM-equation 



1.12) $+(n, .) := ]T K+(n, l)F+(l, .) e t\{Z) 
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Then by use of l|7.2(l and Lemma EO 

^2 $+(n,m)il> q ,-(w,m) = ( ^2K + (n,l)F + (l,m) jip q ,-(w,m) 



V ( — / R + (w)i> + {w,n)^ q . + (w,m)duj(<w) )ip q -(w,m) 
2J {i> q ,-(w, m), R+(w)tfj + (w, n))ip q -(w, m) 



(8.13) = R + (w)ip + (w,n). 



On the other hand, inserting the GLM-equation yields for \w\ — 1 
^2 *+(«, m )'0< ? ,-{w,m) = 



met 



It — L CXJ 

m)i]} q -{w, to) + ^("-, m)K+ 1 (n, n) — K+(n, to) 

m— — oo m=n 
oo g 

- E 7 +)J -^ g , + (pj, O-0«,+(Pi,m) ip q ,-(w, to) 

n-1 

3>+(n, m)ip qt -(w, m) + ip qt -(w, n)K^ 1 {n, n) — ip+{w, n) 



m= — oo 
1 



(8.14) 



j — 1 m— n 



(recall the definition of t/> g ,± from (JfTSJ) and therefore 



(8.15) T(w)h-(w,ri) — ip+(w,n) + R+(w)ip+(w,n), \w\ = 1 

where 

n— 1 



h-(w,n) — 



1 



T(w) \if_|_(n, n) 



E *+(r 



ip q -(w,m) 
ip g -(w,n) 



(8.16) 



j=i 



, W^n-lC^+Cpj),^ ,-(«?)) 



^,-K n )(^H - ft) 

since Green's formula ( [SSI ^ e q- (1-20)) implies for A S cr(H q ) 

oo 

( A -Pi) E i>1,+ (.P3i m )i'<l >-( A > TO ) = -Wn-l(^ g ,+ (p 3 ),^g,-(A)). 
Similarly, we obtain 



h + {w, n) 



ip q>+ (w,n) 



T(w) \K-(n,n) 



E *-(' 



ip qt+ (w,m) 



(8.17) 



E 7 -^-^'' n ) 



i=i 



W n {ij) q -{pj),il) q , + {w)) 

^q,+ ( W , n )(H W ) ~ Pj) 
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with 

n 

$_(n,m)= V" K_(n,l)F~(l,m). 



i=-c 



For n € Z, = 1, we see that /izp(w _1 ,n) = h^(w, n), since K±(n,m) and 
$±(n, m) are real. The functions h^{w,n) are continuous for |w| = 1, w ^ 
w(Ej), since T _1 (iu) is continuous on this set by the Poisson-Jensen formula (|6.25|l 
(|i?±(ui)| < 1 for to ^ w (Ej) by H I7.6I fill and ifj q , T (w 1 m) are continuous on 
dW\{w(nk)}- The functions h T (w, n) have a meromorphic continuation to W\{0} 
with the only possible poles at w(pj) and w([ij). At w(pj) there are no poles, due to 
the zeros of T^ 1 (w) at w(pj). For w — w(pj) we have the same type of singularity 
as ipq,±- in summary, h±(w,n) have simple poles at w{pj) and are continuous at 
the boundary except possibly at w(Ej). 

To study the behavior of h±(w, n) as w — > 0, we recall z^ 1 = —w/a (1 + 0(w)). 
Then 

(~ + 0(v?))W n -i$ q ,+ (p j ),il> q ,-(w)) 

= ( > f W~" +1 ^ q , + (Pj,n-l) + 0(w)), 

(^ + o( w 2 ))^„(^^( P ,),^,+H) 

(— 1)™ f] ._ a g (j) 
and property (B4) implies 

(8.18) V] $±(n,m)iJ qtT (w,m)'ip~^ : (w ) n) = 0(w), w — > 0. 
We conclude that 

(8.19) hm h^(w, n)t/j qt ±(w, n) 



w^o +v ' yY ^'= v ' y T(0)Jf±(n,n) 
H I7.6I (iv) and jOJ imply the following behavior of h T (X,n) as A — > pj 

bm /i T (A,n) = ±7±, i ^±(p i ,n) hm -f yfTTT 1 

2 3 +l 

(8.20) = 7±J V; ± (p j ,ri)(Res p ,T(A))" 1 J] Jp~E h 

1=0 

where h± are defined as in l|6.8|l . 



By virtue of the consistency condition T(w)R + (w) = —T(w)R-(w) we obtain 



h± (w, n) + R± (w)h± (w, n) = 

1 (W) K ' 1 \ w ) v 7 



i/j T (w,n)T(w), \w\ = 1. 
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If we eliminate R T (w) from the last equation and H8.15(l we see 

r 1/2 ( 



T ( w )R2g+2( w )(^+( w i n )' l l>~( w i n ) ~ h+(w,n)h-(w,n)) 



21) 1/2 {h±{w,n)ip±(w,n) - ip±{w,n)h±(w,n)) =: G(w,n), 



^2g+ 2 H 



for \w\ — 1. Observe that G(w,n) — G(w,n) — G(w,n), \w\ — 1, since h±tl)± — 

~ A —1/2 

ip±h± and R2 g +2( w ) are °°-d functions for |w| = 1. The function G(w,n) can be 
continued analytically on W since the difference ip+ip— — h + h- vanishes at the poles 
w(pj) of T(w) by 18.2U[) . Note that the product and hence also h + h_ do not 

have poles at w([ij). Moreover, since W is just the image of the upper sheet, we can 
extend it to a compact Riemann surface W by adding the image of the lower sheet. 
Now by G(w, n) = G(w, n) we can extend G to W by setting G(w, n) = G(w^ 1 , n) 
for \w\ > 1. 

Now let us investigate the behavior at the band edges: If wi ^ w(fij), we obtain 
by H7.30fl . (|8.15|l . and real-valuedness of ip± at the band edges that 



lim i?2g+2HTT( A H ~ A^)M W > n )V'TK™) 

} — *Wi J A A 

3 

= lim <l2U f- H \ ^ + R ± ^ 

= lim Rl ^ +2 U f ~ H) (> ± + 1)* ± + ~ * ± )% 



If wi — w(fij), the same calculation shows that 

lim Rll 2 +2 ( W )T\(K W ) ~ fJ*j)h±(w,n)il>±(w,n) 



= (-l) i+1 C + (n)C_(n) lim <+ 2 H ^\ 1 =0 



It!— »ltJZ 



T(u>) 



by [j73Ujl . where we used -0±(^,n) = i ; C±(n)(A(w) - Mj)~ 1/2 + 0(1). 

Consequently R2 g +2{w)G{w,n) is continuous at ro = wi and vanishes at the 

1 /2 

band edges. Thus the singularities of R 2g+ 2(w)G(w,n) at wi are removable. Fur- 

1 12 

thermore, R 2g+ 2(w)G(w,n) is purely imaginary for \w\ = 1 and real on the slits 
and hence must vanish at wi by continuity. So the singularities of G(w, n) at wi 
are removable as well. Thus G is holomorphic on all of W and vanishes at w = 0, 
that is, G(w,n) = which implies (compare (B4)) 

lim ( ip+(w, n)ip^ (w, n) — h+(w, n)h- (u>, n) J 
= K + (n,n)K-(n,n) - (T(0) 2 K + (n, n)K^(n, n))' 1 = 0. 

Using |jH2Il we finally obtain from T(0) 2 = (K + {n,n)K-(n,n)y 2 that 
(8.22) a+{n) = a-(n) = a(n), Vn E Z. 
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It remains to prove b+(n) = b-(n). Proceeding as for G(w, n) we can show that 
T(w)R 2 gi^(w)( K '4!j r (w, n)-0_(w, n + 1) — h + (w,n + X)h-(w, n)) = 

UAKw) - m) , 

(8.23) — ~~m (h + (w,n+ \)ip + (w,n) - ip + (to, n)h+ (w, n + 1)) 

R 2g+2( W ) 

is a constant equal to — l/o(n). Thus 

W(w,n) := a(n) (ip+(w, n)'ij)—(w,n + 1) — h+(w,n + (w,n)) 

(8 94) = R 2 9 l2M 

( ' ] rHn 3 (AH- ft )- 

and computing the asymptotics at w = (compare (|4.3|) ) we see 
(8.25) = W{w,n) - W{w,n- 1) = A+(0)A_(0)(6+(n) -6_(n)) 

and in particular = b_(n) = b(n). 

Our operator has the correct norming constants since as in (|6.12|l it follows 

2g+i 

.26) ^2i> + ( Pj ,n)^( Pj ,n) = (Res Pj T(A))" 1 ]^[ VpT^ 



r 1 



and by 

^± (ft ' n )^± (Pi ' n ) = T±j ■ 

n£Z 

□ 
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